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ABSTRACT 

An  analysis  and  test  rig  measurements  of  the  per- 
formance of  a  transonic  axial  turbine  are  reported.   The 
purpose  was  to  confirm  the  accuracy  of  measurements  made 
in  a  test  rig  which  was  designed  to  separate  the  losses 
occurring  in  the  stator  from  the  losses  occurring  in  the 
rotor  blade  rows.   The  analysis  was  programmed  for  the 
Hewlett-Packard  21-MX  computer.   Reasonable  agreement 
between  predicted  and  measured  characteristics  was  obtained 
using  experimentally  determined  losses  in  the  computer 
program.   Lack  of  agreement  was  noted  using  theoretical 
values.   It  was  concluded  that  the  rotor  was  not  choked 
at  the  conditions  in  the  tests,  and  that  the  test  rig 
measurements  were  valid.   A  successful  technique  for 
smoothing  the  data  obtained  from  the  rig  is  also  reported. 
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I.   INTRODUCTION 

The  transonic  turbine  test  rig  installation  at  the 
Naval  Postgraduate  School  Turbopropulsion  Laboratory  was 
designed  to  study  the  effects  on  turbine  performance  of 
varying  axial  and  tip  clearances,  to  study  the  effects  of 
blading  design  on  turbine  performance,  and  to  allow  the 
separate  determination  of  stator  and  rotor  losses  in  an 
operating  machine . 

Until  the  work  of  Solms  (Ref.  1)  the  separation  of 
rotor  and  stator  losses  through  test  rig  measurements  had 
not  been  attained  satisfactorily.   Through  improvements 
in  hardware  and  instrumentation  and  improvements  in  the  data 
reduction  process  the  stator  and  rotor  losses  were  deter- 
mined separately  and  were  reported  in  Ref.  1. 

Anomalies  remained,  however.   Specifically,  the  turbine 
coBif figuration  designated  as  Turbine  C  in  Ref .  1  gave 
different  results  when  compared  in  terms  of  "referred" 
quantities  depending  on  whether  the  discharge  was  to 
atmospheric  pressure  or  to  a  region  of  reduced  pressure. 
(Turbine  C  had  converging-diverging  stator  passages  in  an 
axial  entry,  single  impulse  stage.   The  turbine  was  designed 
to  operate  in  the  transonic  range.)   In  addition,  R&f.  1 
reported  considerable  scatter  in  the  loss  coefficients. 

In  the  work  of  Robbins  (Ref.  2)  it  was  shown  that 
discharge  pressure  affected  the  measurements  of  flow  rate 
into  the  stage  and  also  affected  the  labyrinth  leak  rate. 


Accordingly,  Robbins  determined  accurately  the  flow  rate 
into  the  stage  and  the  leak  rate  through  the  labyrinth  seal 
(See  Fig.  1)  for  all  operating  conditions.   The  results 
then  obtained  for  turbine  C  were  reported  fully  in  Ref .  2. 
The  continued  presence  of  scatter  in  the  measured  loss 
coefficients  was  reported  and  a  smoothing  technique  to 
eliminate  the  scatter  was  suggested. 

Before  the  Turbine  Test  Rig  could  be  used  to  measure 
the  effects  of  varying  parameters i 

1.  The  scatter  in  the  loss  coefficients  had  to  be 
eliminated. 

2.  The  overall  accuracy  of  the  performance  results 
evaluated  from  the  rig  measurements  had  to  be 
verified  in  some  way. 

The  resolution  of  these  problems  was  the  goal  of  the 
present  work  and  is  the  subject  of  this  report.   First, 
a  satisfactory  method  was  found  for  smoothing  the  loss 
coefficients.   The  method  is  described  in  Section  III. 

The  approach  taken  to  verify  the  performance  of  the 
rig  was  to  first  devise  an  analysis  which  predicted  the 
performance  of  the  turbine  in  terms  of  unknown  loss  coeffic- 
ients, and  then  to  show  that  the  measured  loss  coefficients 
were  consistent  with  the  predicted  behaviour. 

A  description  of  the  Turbine  Test  Rig  is  given  in 
Section  II.   The  analysis  of  the  behaviour  of  the  turbine, 
involving  a  comparison  of  an  analytical  prediction  with 
the  results  of  a  short  test  program,  is  described  in 
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Section  IV.   Details  of  the  analysis  and  the  computer 
program  are  given  in  Appendix  A. 
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II.   TURBINE  TEST  RIG  INSTALLATION 

A.    DESCRIPTION 

The  test  installation  consists  of  three  major  compon- 
ents? an  Allis-Chalmers  twelve  stage  axial  flow  compressor, 
an  exhauster  assembly,  and  the  turbine  test  rig  (TTR) 
itself. 

The  compressor  is  the  source  of  driving  air  for  the 
TTR  and  for  the  exhauster  assembly.   Fig.  2  shows  the 
piping  arrangement.   Turbine  air  passes  through  the  first 
settling  tank  into  an  eight-inch  pipe  containing  a  flow 
nozzle,  into  the  second  settling  tank  and  into  the  turbine. 

Fig.  3  shows  the  plenum,  the  floating  stator  assembly, 
the  rotor,  and  the  dynamometer  (Ref.  1).   Pressure  ratios 
of  6*1  can  be  achieved  when  the  system  is  hooded.   The 
hood  was  needed  to  achieve  high  pressure  ratios  in  the 
tests  reported  here.   Fig.  k   shows  the  turbine  blading 
of  the  stator  and  rotor.   Ref.  3  contains  detailed  descrip- 
tions of  the  test  rig  hardware. 

The  floating  stator  assembly  shown  in  Fig.  3  permits 
measurements  of  the  axial  force  and  the  torque  on  the 
assembly.   Axial  and  rotational  movements  are  constrained 
by  calibrated  force  transducers  that  are  heat  insensitive. 
These  measurements,  together  with  wall  static  pressure 
measurements,  allow  the  determination  of  the  average 
axial  and  tangential  velocity  components  at  the  stator  exit. 

In  this  report  one  configuration  designated  Turbine  C 
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was  tested,  the  geometry  of  which  is  shown  in  Fig.  5. 
Table  I  describes  the  geometry  quantitatively.   The  stator 
blade  profile  is  shown  in  Fig.  6.   The  blades  of  the  stator 
generate  a  converging-diverging  nozzle  shape.   Pressure 
measurements  were  taken  at  the  locations  shown  in  Fig.  6. 
The  pressures  necessary  to  the  analysis  of  the  stator  axial 
force  were  taken  at  the  locations  shown  in  Fig.  7. 

B.  TEST  MEASUREMENTS  AND  ACCURACY 

1 .  Mass  Flow  Rates 

Appendix  A  of  Ref .  2  gives  a  detailed  description 
of  the  method  used  to  determine  both  the  turbine  flow  rate 
and  the  labyrinth  seal  leak  rate. 

2.  Forces.  Torques.  Temperatures,  and  Pressure 
Ref.  3  and  Ref.  5  give  calibration  procedures 

for  the  TTR.   Identical  procedures  were  employed  here. 
Table  II  of  Ref.  1  gives  the  expected  accuracies  of  the 
measurements . 

C.  TESTING  AND  DATA  REDUCTION 

The  TTR  data  collection  system  is  described  in  Ref.  k. 
Appendix  D  of  Ref.  1  gives  a  detailed  explanation  of  the 
turbine  test  procedures.  Those  procedures  were  followed 
here  with  the  exception  that  a  constant  RPM  was  held  and 
the  pressure  ratio  varied  over  the  desired  range.  The 
data  reduction  method  developed  in  Ref.  1  and  Ref.  2  was 
revised  as  described  in  Appendix  B. 
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Ill,   DATA  SMOOTHING  TECHNIQUE 

The  sensitivity  of  the  loss  coefficients  to  variations 
in  measured  quantities  was  shown  in  Ref.  2.   It  was  also 
stated  that  some  variation  is  unavoidable  since  measure- 
ments are  taken  over  a  period  of  more  than  a  minute.   It 
was  also  pointed  out  in  Ref.  2  that  the  parameter  most 
important  to  the  calculation  of  the  loss  coefficients  was 
P.  (the  average  pressure  at  the  stator  exit).   P.  is  an 
average  pressure  that  can  not  be  measured  directly.   It 
is  derived  as  described  in  Ref.  1  from  many  other  measure- 
ments. Significant  scatter  was  observed  in  the  variation 
of  P.  as  speed  was  varied  at  fixed  pressure  ratio.   However, 
it  was  found  that  the  hub  and  tip  pressures  (P.  and  P.) 
measured  just  downstream  of  the  stator  varied  smoothly  over 
the  same  range.   Since  these  two  pressure  were  measured 
directly,  and  since  they  varied  smoothly,  it  was  assumed 
that  the  pressure  behind  the  stator  must  vary  smoothly  also. 
Consequently,  it  was  determined  that  a  polynomial  curve  fit 
could  be  used  to  describe  the  variation  of  P.  as  speed  was 
varied.   The  variation  of  P.  was  represented  as  a  function 
of  P.  ,  P.,  and  the  isentropic  head  coefficient  (K.  ) 
in  the  formi 

°~  =  iVf^  =  /1-  +A>Kit+AiKixs-'-.  (i) 

Where  A^  is  the  polynomial  coefficient. 

Then,  P1/^to  (where  P^  is  total  pressure  upstream  of 
the  stator)  was  computed  using  the  expression 
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pjl  =<rtoi)+  £JL 


P  \     P  /IP 

*to      \    to  /    \  toy 

Fig.  8  is  an  illustration  of  this  procedure,  which 
was  added  to  the  "bulk  process"  data  reduction  program. 
The  data  points  in  Fig.  8  were  taken  from  Runs  6  and  7 
in  Ref.  2.   The  two  Runs  were  for  the  same  conditions  but 
were  made  at  different  times.   It  can  be  seen  that  the 
trends  are  the  same  for  both  runs  but  the  scatter  is 
considerable.   The  scatter  is  the  cause  of  the  scatter  in  the 
calculated  loss  coefficients.   The  lines  on  Fig.  8  are  the 
polynomial  approximations  according  to  Eq.  (1)  for  the  two 
runs.   It  can  be  seen  that  the  polynomial  approximation 

averages  the  data  and  maintains  the  original  trend.   The 

st 
chosen  smoothing  function  was  a  1   degree  polynomial. 

Fig.  9  shows  the  stator  loss  coefficient  prior  to 
smoothing  for  the  points  in  Run  #7.   Fig.  10  shows  the 
same  data  after  smoothing. 

Similarly,  Fig.  11  shows  the  rotor  loss  coefficients 
for  Run  #7  prior  to  smoothing  and  Fig.  12  shows  the  same 
data  after  smoothing. 

It  can  be  seen  from  these  figures  that  the  scatter 
has  been  removed.   This  is  particularly  true  for  the  rotor 
loss  coefficients.  As  pointed  out  above,  it  is  believed 
that  the  observed  scatter  was  due  to  the  sensitivity  of 
the  loss  coefficients  to  small  changes  in  measured  quan- 
tities during  the  data  collection  process.   The  smoothing 
technique  removes  the  random  variations  recorded  during 
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data  collection  and  results  in  a  much  more  realistic 
representation  of  the  variation  in  the  losses. 

The  smoothing  technique  was  incorporated  into  the 
data  reduction  program  as  described  in  Appendix  B. 
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IV.   ANALYSIS  OF  TURBINE  PERFORMANCE 

A.  APPROACH 

In  order  to  determine  if  the  performance  evaluated 
from  the  rig  measurements  was  accurate  an  analysis  to 
predict  the  behaviour  of  the  test  turbine  was  carried  out 
and  programmed  in  BASIC  language. 

A  performance  test  was  then  conducted  in  a  particular 
way  in  order  to  provide  a  comparison  of  the  measured  with 
the  predicted  behaviour. 

B.  ANALYTICAL  PREDICTIONS 

The  flow  through  the  turbine  was  analysed  using  a 
pseudo-1  Dimensional  compressible  approach.   The  analysis 
is  described  in  detail  in  Appendix  A,  together  with  the 
computer  program  which  was  used  to  obtain  predictions  of 
the  turbine  performance.   One  of  the  inputs  which  the 
program  requires  is  the  rotor  passage  loss  coefficient  at 
zero  incidence.   Using  the  method  given  by  Vavra  in  Ref.  6 
it  was  determined  that  the  rotor  loss  coefficient  should 
have  the  value  .251^.   However,  results  of  previous  tests 
of  the  turbine  indicated  that  the  rotor  loss  coefficient 
was  rarely  as  high  as  .251^  and  could  be  as  low  as  .1. 
Therefore,  the  performance  of  the  turbine  was  analysed 
using  rotor  loss  coefficients  of  .251^  and  .1.   The  pre- 
diction program  was  run  for  both  15,000  and  18,000  RPM  with 
an  assumed  rotor  loss  coefficient  of  .251^,  and  for  18,000 
RPM  with  an  assumed  rotor  loss  coefficient  of  .1 
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The  above  parameters  were  chosen  to  obtain  a  prediction 
of  the  turbine  performance  in  a  range  in  which  experimental 
data  could  be  obtained.   In  particular,  the  analysis  could 
be  used  to  predict  the  pressure  ratios  at  which  choking 
would  occur  in  the  rotor  as  well  as  in  the  stator.   The 
pressure  ratio  at  which  choking  occurred  in  the  rotor  could 
then  be  established  experimentally  in  a  test  conducted  at 
fixed  speed.   An  examination  of  the  choking  condition  was 
considered  to  be  a  first  test  of  the  performance  analysis. 

The  results  of  the  analysis  for  the  parmameters  given 
above  are  shown  in  Fig.  13,  Fig.  Ik,    and  Fig.  15.   What  is 
shown  in  the  figures  is  the  map  of  the  range  of  values 
which  unknown  parameters  in  the  analysis  can  have,  that 
leads  to  a  solution. 

Figures  16-19  show  predicted  performance  parameters 
for  the  case  of  18,000  RPM  and  assumed  rotor  loss  of  .1. 
These  results  will  be  discussed  in  conjunction  with  the 
results  of  the  turbine  test  run. 

C.    EXPERIMENTAL  TURBINE  TEST 

In  order  to  examine  the  occurrence  of  rotor  choking 
the  turbine  was  run  at  constant  RPM  and  the  pressure  ratio 
across  the  stage  was  increased  in  increments  by  lowering 
the  back  pressure.   The  point  at  which  the  horsepower  ceased 
to  increase  for  an  increase  in  pressure  ratio  was  examined 
to  determine  the  choking  point.   At  the  condition  where 
the  flow  reaches  a  Mach  number  of  unity  at  the  exit  of  the 
rotor,  the  power  produced  by  the  turbine  can  not  be  changed 
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by  altering  the    downstream  pressure. 

Tests   were    conducted  in  this   manner   at   15 f 000  and 
18,000  RPM.      The   controlled  parameters   for   the   tests   are 
given  in  Table   II.      The   reduced  data  is  given  in  Table    III. 
The   referred  horsepower   is   shown  plotted  versus   the   pressure 
ratio   for   the    15,000  RPM  run  in  Fig.    20  and  for   the   18,000 
RPM  run  in  Fig.    21. 

The   results   are    discussed  in  the   next  section. 
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D.    COMPARISON  OF  ACTUAL  AND  PREDICTED  PERFORMANCE 

Figures  13 1  l^+i  and  15  are  maps  of  possible  solutions 
for  the  flow  through  the  turbine,  when  the  parameters  which 
are  unknown  are  allowed  to  vary.   It  is  noted  first  that 
the  predicted  range  of  solutions  extends  to  pressure 
ratios  (Py^to^  Del°w  *^e  predicted  choking  line.   The 
explanation  for  this  apparent  anomaly  is  that  the  program 
calculates  all  possible  solutions,  and  for  any  point  below 
the  line  a  particular  combination  of  losses  and  blockage 
factors  existed  which  would  allow  a  solution  at  that  point 
without  choking.   It  should  be  pointed  out  that  any  point  on 
the  plot  can  be  brought  to  the  choking  line  by  reducing 
the  stator  loss  by  a  very  small  amount.   This  is  the 
procedure  that  was  followed  to  get  the  range  of  values 
at  choking  shown  in  Figures  16-19. 

Fig.  13  is  the  map  produced  by  the  program  for  RPM 
equal  15»000  and  for  an  assumed  rotor  passage  loss  coeffic- 
ient equal  to  .251**.   As  can  be  seen  from  the  figure  the 
pressure  ratio,  (Pp^to^*  at  choking  was  about;  .28,  cor- 
responding to  a  stage  pressure  ratio  (P+q/Po)  °^  3*57. 
Fig.  Ik   is  the  map  for  RPM=18,000.   Note  that  the  predicted 
stage  pressure  ratio  for  choking  is  again  about  3«57» 

Fig.  20  and  Fig.  21  show  the  variation  of  the  referred 
horsepower  vs.  the  pressure  ratio  which  were  measured  in 
turbine  tests  at  15,000  and  18,000  RPM  respectively.   It 
can  be  seen  that  the  referred  horsepower  did  not  become 
independent  of  the  back  pressure  at  the  pressure  ratios 
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achieved  in  these  tests.   It  was  concluded  that  the  stage 
did  not  choke  at  the  predicted  pressure  ratio  of  3*57* 
It  was  noted  however  that  at  the  highest  pressure  ratios 
obtained  at  18,000  RPM,  the  slope  of  the  horsepower  curve 
was  becoming  smaller. 

The  velocity  diagrams  in  Figures  23-31  are  also  consis- 
tent with  the  argument  that  the  rotor  did  not  choke  during 
the  turbine  tests.   The  velocity  from  the  rotor,  V"2, 
behaves  smoothly  and  in  a  predictable  manner  throughout 
the  pressure  range  at  both  test  speeds.   This  might  not 
be  expected  if,  following  choking,  shock  waves  appeared 
downstream  of  the  rotor  exit  plane. 

Note  the  values  of  rotor  loss  coefficient  given  for 
the  test  results  in  Table  III.   With  the  exception  of 
point  #10,  which  was  considered  to  be  in  error,  they  were 
all  smaller  than  the  value  (.251*0  assumed  in  the  first 
performance  calculations.   As  explained  in  Appendix  A, 
the  value  of  the  rotor  passage  loss  coefficient  used  in 
the  program  is  the  smallest  value  that  can  be  calculated 
for  the  overall  rotor  loss  coefficient.   Therefore,  if 
the  experimental  results  were  accurate,  the  computer  program 
could  not  predict  the  correct  performance  since  it  could 
never  calculate  a  rotor  loss  less  than  the  input  value  of 
the  passage  loss,  which  was  .251^.   It  was  as  a  consequence 
of  this  observation  that  the  performance  was  re-calculated 
using  a  rotor  passage  loss  coefficient  equal  to  .1. 

Figures  16-19  show  the  choking  behaviour  of  the 
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turbine  predicted  using  the  computer  program  for  an  assumed 
rotor  passage  loss  of  .1.   In  these  figures,  the  data  for 
test  point  #6  is  also  shown.   In  this  case  it  can  be  seen 
that  the  predicted  pressure  ratio  (Pp/^to^  a^  choking  is 
about  .2^+5»  corresponding  to  a  stage  pressure  ratio  (P+0/P2^ 
of  about  ^.08.   The  latter  is  slightly  higher  than  the 
maximum  pressure  ratio  that  was  attained  in  the  turbine 
tests.   (At  18,000  RPM  the  highest  attainable  pressure 
ratio  was  3«97»  due  to  the  characteristics  of  the  dynamom- 
eter. ) 

Fig.  1?  is  the  range  of  efficiencies  predicted  for  the 
choked  condition.   It  is  noted  that  the  efficiency  measured 
at  test  point  #6  was  reasonably  close  to  the  predicted 
maximum  efficiency. 

Fig.  18  shows  the  range  of  stator  losses  for  which 
solutions  existed,  in  comparison  with  the  value  of  stator 
loss  measured  at  test  point  #6.   It  can  be  seen  that  the 
measured  value  intersects  the  predicted  range  of  possible 
solutions.   Fig.  19  shows  the  predicted  range  of  the  rotor 
coefficient.   Again,  it  can  be  seen  that  the  measured  value 
at  test  point  #6  overlaps  the  predicted  range  of  possible 
solutions. 

Fig.  16  shows  the  predicted  range  of  horsepower 
compared  with  the  horsepower  measured  at  test  point  #6. 
The  measured  horsepower  was  slightly  greater  than  the  max- 
imum value  which  was  predicted. 
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E.    DISCUSSION 

The  results  shown  in  Figures  16-19  illustrate  the 
uncertainty  in  the  prediction  of  the  performance  of  the 
turbine  when  only  the  rotor  passage  loss  and  blockage 
factor  are  known.   It  is  recalled  that  the  analysis 
satisfies  only  continuity  through  the  stage,  and  values  of 
an  additional  loss  coefficient  and  an  additional  blockage 
factor  must  be  established  to  obtain  a  unique  solution. 
It  is  the  aerodynamic  shaping  of  the  surfaces  which  deter- 
mines these  factors.   Figures  16-19  show  the  possible 
range  of  performance  that  results  simply  from  the  areas 
of  the  passages  and  the  blade  angles. 

The  irregular  shapes  of  the  bounds  on  the  possible 
solutions  in  Figures  16-19  are  interesting.   There  is  no 
obvious  explanation  for  the  reduced  range  of  solutions  near 
kbl  =  .8!. 

The  results  obtained  in  the  present  work  suggest 
that  the  performance  of  the  turbine  as  measured  in  the 
turbine  test  rig  is  reliable.   It  had  been  thought  prev- 
iously that  the  rotor  loss  measurements  were  too  low. 
Here,  an  analysis  was  carried  out  to  predict  the  performance 
of  the  turbine,  and  the  predicted  results  have  shown  very 
good  agreement  with  the  performance  measurements  in  the  test 
rig. 

Whether  or  not  the  rotor  was  choking  had  also  been  a 
question  in  the  past.   In  the  present  work  both  the  analysis 
and  experimental  tests  have  shown  that  the  rotor  was  not 
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choked  in  any  test  to  date.   The  computer  program  with  an 
assumed  rotor  loss  of  .1  has  shown  that  the  test  rig  should 
not  be  choked  at  any  pressure  ratio  below  about  ^.08.   Test 
data  has  shown  that  the  rotor  was  not  choked  at  a  pressure 
ratio  just  slightly  below  ^.0. 

It  is  also  of  interest  to  note  that  for  a  rotor 
passage  loss  coefficient  of  .251^  the  predicted  pressure 
ratio  was  shown  to  be  independent  of  RPM.   At  both  15,000 
and  18,000  RPM  the  choking  pressure  ratio  (P.  /P?)  was 
calculated  to  be  about  3«57-   More  results  are  needed  to 
confirm  that  the  choking  pressure  ratio  is  indeed  indepen- 
dent of  speed. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

1.  The  rotor  of  Turbine  C  was  not  choked  at  any  pressure 
ratio  tested  so  far.   This  conclusion  is  based  on  the 
prediction  of  the  computer  program  and  also  on  the 
turbine  test  results. 

2.  The  loss  coefficients  measured  in  the  TTR  and  smoothed 
as  described  in  this  report  can  be  aocepted  as  being 
truly  representative  of  the  losses  in  the  separate 
blade  rows  of  the  turbine . 

3.  In  particular,  the  results  from  the  computer  program 
suggest  that  the  magnitude  of  the  measured  rotor  loss 
coefficients  is  probably  correct.   The  method  used 

by  Vavra  in  Ref .  6  predicts  much  larger  values  of  loss 
coefficients  for  the  present  rotor  geometry  than  those 
which  were  measured.   When  the  measured  value  of  the 
rotor  loss  coefficient  was  entered  into  the  computer 
•    program,  there  was  good  agreement  between  all  the 
calculated  and  the  measured  turbine  performance 
parameters.   When  the  higher  (calculated)  value  of  the 
rotor  loss  was  entered,  there  was  a  pronounced  dis- 
agreement between  the  calculated  and  the  measured 
performance. 

k.        The  computer  program  should  be  used  (and  progressively 
developed)  in  conjunction  with  all  tests  carried  out 
in  the  turbine  test  rig.   Most  importantly,  a  test 
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should  be  conducted  to  determine  the  rotor  choking 
condition  experimentally.   This  will  require  the  purchase 
of  a  water-brake  dynamometer  to  extend  the  power-speed 
range  of  the  test  rig. 

5.   Since  the  accuracy  of  the  test  rig  measurements  is 

no  longer  in  doubt,  experiments  to  determine  the  effect 
of  parameter  changes  (e.g.  axial  and  tip  clearances) 
on  turbine  performance  can  go  ahead. 
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TABLE  I 
TURBINE  GEOMETRY 


TURBINE 

STATOR 

ROTOR 

C 

1 

1 

DESCRIPTION 

A*   2 
(in? 

NUMBER 
OP 
BLADES 

*M 

(in.) 

STATOR  1 

CONVERGING 
DIVERGING 
NOZZLES 

2.9058 

31 

4.184 

0.5775 

ROTOR  1 

CIRCULAR 
ARC 

7.119 

60 

Rm1=4.l93 
Rm2=4.250 

h^O.732 
h2=0.8475 
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TABLE  II 
PARAMETERS  FOR  TURBINE  TEST 


POINT  # 

RPM 

Pto/P2 

1 

15,000 

2.01 

2 

15,000 

2.^2 

3 

15,000 

2.93 

k 

15,000 

3.^9 

5 

15,000 

3.68 

6 

18,000 

3.9^ 

7 

18,000 

3.^+7 

8 

18,000 

2.98 

9 

18,000 

2.52 

10 

18,000 

2.02 
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TABLE  III  TURBINE  TEST  RESULTS 


POINT  #1 


PRESSURE  RE TIG 

REF  FLOW  RATE 

REF  RPM 

REF  ROTOR  MOMENT 

REF  HP 

EhF.  T-S 


2 .  0 1 1 3  5  3  3  3  S 
1 .  8 £889 0888 
14262. 15617 
102,,  5192926 

•"' . "!  1  4  '-^ ! "! '. "!  1  r*!  '■**' 

0. 712723663 


STAT OR  LOSS  THEOR 
STAT OR  L033  COEFF 
ROTOR  LOSS  COEFF. 
ROTOR  LOSS  THEOR. 

I  SEN.  HERD  COEFF,. 


-  8„ 137100795 

=  0.220295023 

=  0. 193137353 

=  0.  3  £6  131836 


4  „  0362^0 ,-" 


T  H .  I'  E  G .  0  F  R  E  A  C  T 1 0  N  =  0.825  3  8  664  4 

A  C  T  .  I)  F  i_, .  i ;:  F  F:  E  A  C  T I  0  N  =  -  8 .,  O  8  4  6  7  3  1  8  2 


VELOC 

I  TV 

TRIANGLE 

DATA 

V  1  :~ 

973 

.  f  •„'  ■-'  iC  ( 

V  2  = 

c.  U  4  ■  ii  3  -j  -'  -t 

Vfll  = 

2  26 

, 65467 

VA2  = 

1  .j  -j  .  4  c  £.  a  1 

y  1 1  •!  -■ 

946 

. 93633 

VU2= 

1 3  2 .  4  3  9  3  2 

fl-LPHf 

1  1=' 

76. 533 

37    AL- 

p h  R  2  =    :i  8  . 

4473 

BETH 

■I 

68. 444 

IO    BE 

TA  2  =   -69 „ 

3457 

POINT  #2 


PRESSURE  RATIO 

REF  FLOW  RATE 

REF  RPM 

REF  ROTOR  MOMENT 

REF  HP 

EFF.  T-S 


=  2.42 3  5 2    ~ 8 s 

-  1.0  &  8  0  0  0  8  0  0 
=  1  4  3  3  "■' .  O  5 ::'[  3 
=  132  ,,334797 

~  SO ,,04433  335 

-  0 .,  7  4  3  9  7  3  3  5  8 


ST AT OR  LOSS  THEOR 
ST AT OR  LOSS  COEFF 
ROTOR  LOSS  COEFF,, 
ROTOR  LOSS  THEOR. 


I  SEN,.  HEAD  COE 
TH.  DEC, 
ACT.  DEC. 


•-  0.  1O7100795 

=  0. 133632754 

=  0. 133537254 

=  0.455631440 


=  4,941765312 
REACTION  =-8. 76562 E- 03 

F  R  E  A  C  T 1 0  N  =  -  8 .  2  7  2  7  9  6  9  0  1 


VELOCITY  TRIANGLE  DATA 


VI  = 
VfiJ> 

vui- 

1126 

c.  b  '0 
1.094 

06734 

4  €  -1 4  9    ' ,' 

32844 

2  =   134.36934 
A 2=   I73,i  93521 
U2=    44.23066 

ALPHl" 
BETA 

1 = 
•i 

76  .  ,:■  6  4 '■■::'.. 

6  4  .  0  2  3 1 Z 

ALPHA  2=    13 
BETA  2  =   -78 

„  7  9   -1  6  i:: 'I 
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TABLE  III  TURBINE  TEST  RESULTS  (CONT. ) 
POINT  #3 


P  R  E  3  S  U  R  E  R  A  T  1 0  =  2 .  9  3  2  5  7  8  7  4 

F:  E  F  F  L  0  N  F:  R  T  E  = .  1 .  9  2  8  @  8  0  0  0  8 

REF  RPM  =  14297.66424 

REF  ROTOR  MOMENT  =  161.7256537 

REF  H  P  ~-    3  6  u  6  y  6  8  5  4  7 

EFF„  T-S  =  0.769522864 

S TRIOR  LOSS  THEOR.  =  '0. 107180795 

S  T  fl T 0  R  L  0  S S  C 0  E  F  F .  =  0.1 7  2  8  6  5  7  8 7 

F  0  T  i j  R  L  0  a  3  C  0  E  F  F  .  =  8 „  0  3  3  5  2  f  ■  2  2  3 

F:  0  T  0  R  L  i  J  8  S  T  H  E  0  R ,,  =  8 .  4  3 179246 

= -  @ .  8  3 1  2  3  9 1 8  9 
T 1 0 N  =  -0.241795592 

VELOCITY  TRIANGLE  DATA   ' 

VI  =  1248.51552  V2  -  213,93763 
V  A 1  =  3 1  4  .  7  1 8  9  5  V  A  2  =  2  0  6 .  8  8  3  7  5 
V U 1=120 8 . 1 9  8 2  4    V U 2 =   - 7 1 . 6  4 3 5  6 

H  L  !■  H  H  1  -    ?  5  .  3  9  ?  S  3    R  L  F  H  R  2 =   -  i "} ,  3  $  9  9  !■  :i 
B  E I A  1  =    6  4 .  4  9  9  8 1    B  E  "i"  A  2  =   -71.  7  8  6  0  7 


IS 

;EN 

.    HERD 

:oeff 

TI- 

l. 

DEG.    OF 

RE  AC 

T  I  0 

AC 

:T. 

DEG.    0 

r    p£f 

CTI 

POINT  #4 


PRESSURE  RATIO  =  3.490951547 

REF  FLOW  RATE  =  1 . 823888008 

REF  RPM  =  14298.27554 

R  E  F  R  0  T  0  R  M  0  M  E  N  T  =  1 3 1.571 3  7 1 2 

REF  HP  =  41. 16922137 

EFF.  T-S  a  0,,  753539626 

S T  fl  T 0 R  L 0  S S  T H  E  0 R .  =0=1 8 7108  ? 9 5 

S T A T 0  R  L  0  S S  C 0  E  F  F .  =  8.15  8 8  6  8 2  5  9 

R 0  T 0 R  L  0 3  S  C 0  EFF.  =  8.1 3 4  4  9 1 1 8 5 

ROTOR  LOSS  THEOR,  =  0.427987152 

I S  E  N .  H  E  A D  C 0  E  F F .  =6.70 3 6 6 9  5 3  5 
TH.  DEG.  "OF  REACTION  =-8.814634533 
R C T .  HE G .  0 F  R E A C T 1 0 H  =  --3.21 3 8 8 4 7 3 8 

VELOCITY  TRIANGLE  DATA 

' , !  A 1  =   3  5  5  .  9  8  6  3  3    '■ '  A  2  =   2  3  8 .  2  6  8  5  5 
V  U 1 s  1  2 8 5 . 8 1 4  2 8    V  U 2 =  -151. 4 3 6 3 8 

F;  L  F:  !•  I A  1  =    7  4 .  5  5  6  3  1    A  L  P  H  A  2  =   -  3  2 .  ■■!•  3  9  6  9 
B  E  T  A  1  =    6  4 .  2  6  9  1  2    B  E  "i"  A  2  =   -  7  I  .  4 1  3  2  3 
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TABLE  III  TURBINE  TEST  RESULTS  (CONT.) 


POINT  #5 


PRES 

:.:.  Ij 

RE 

RATIO 

REF 

-L 

3W 

-ATE 

REF 

:;:p 

•1 

REF 

?  fj 

TOR 

11  ONE  NT 

REF 

-IP 

EFF. 

T 

-  S 

STfl'T 

JR 

LU 

SS  THE OF 

•   ST  AT 

:ir 

LU 

3S  CO EFF 

ROTO 

'"' 

L..GS 

3  CO EFF. 

ROTOR 

-OS 

1  THE OR. 

J.  bi-  .'J  3  9  b  DO 

1  .  O,;:  000 OOt U 


I  SEN.,  HERD  CO  EFF. 


=  l 


7432 


=  i  S  6  «  3  6  8  2  4  3 6 
■  i-_  ii  .''  i  f  6  '■$  -j  fe 

••-  0.752S03647 

:~  O.  i  07 10O795 
=  0. 140979499 
;--  O.  192867968 

f.|   -1  •"''  C'  A  -"i  •-.  i~i  a  ~? 


6.  9505-49' 


T  H .  D  E  G .  0  F  R  E  fl  C  T 1 0  N  =  8  „  9  O  8  8  2 1  -  8 : 
A  C  T  .  IJ  E  G .  0  F  R  E  A  C  T 1 0  N = -  0 , 1  2  3  €  3  675' 


'El 


I  TV  TRIANGLE  DATA 


VI  =  1362. 18537 
VA1=  361. 18549 
VUi=  1313.45034 


'A2= 


701 ( 42  ?9S 
250. 0 @96 3 
168, 39007 


ALPHA  1= 
BETA  1.  = 


7  4  i,  6  2 1  ■'  5  •_ 
h  4  B  g  9  y  g ;_ 


ALPHA 


hi  2  = 


9  b  J.  b  i 
0212: 


POINT  #6 


PRESSURE  RATIO 

REF  FLOW  RATE 

REF  RPM 

REF  ROTOR  HONE  if 

REF  HP 

EFF.  T-S 


■.:• .  -1  j  o  .j  i  c  4  £  r 

1  .  020000000 

17091 . 82469 
169,  22:.  2167 
45.39160403 
O. 784486361 


ST ATOP  LOSS  THEOR 
STflTOR  LOSS  CO EFF 
ROTOR  LOSS  CO EFF. 
ROTOR  LOSS  THEOR. 

I  SEN,  HEAD  CO EFF. 


=  0.1 0 7 1 0  O  7  9  5 

=  O. 141347539 

;-  8,1 8  7  7  2  2  7  4  0 

=  0. 421689292 

=  5.850  9  O  2  8 ' 


T  H  .  I)  E  G .  0  F  R  E  A  C  T 1 0  N  =  O  „  1  04266156 
fl C T .  BEG.  0  F  R E fl C T 1 0 N =  9.175212 9  3 7 

VELOCITY  TRIANGLE  DATA 


VI  =  13  IS. 03669 
Vfil  =  334, 23319 
VU1=  1274, 94005 

ALPHA    1=         75„3i 
BETA    1    =        61.5: 


V  II  ■-.  - 


.:  t  9 , 

i"4  4  -:_ 

o 

:•  6  0  . 

3475 

:::, 

-  f  9 « 

3348 

.j 

AL 

:'Hfl 

2  = 

..-  ! 

5  .  0  C: 

BE 

FA  2 

:~ 

— ;i 
~°  f 

0 ,  5  i 

i  >•+ .-:  r- 
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TABLE  III  TURBINE  TEST  RESULTS  (CONT. ) 


POINT  #7 


PRES 
REF 

REF 
REF 

REF 
EFF. 


SURE 

FLOW 

RPM 

ROTO 

HP 


:ATIO 
IflTE 

MOMENT 


•;=  3  „  46"'" 553  140 
=  1. .  820088066 
=  17039.50315 
=  155. 9 6  S 0 S 4  2 
=  42. 16691707 
=  0.734673443 


STMT  OR  LOS: 
S TRIOR  LOS: 
ROTOR  LOSS 
ROTOR  LOSS 


COEF 
COEFF 

THEOP 


=  0. 13710O795 
=  3.165054929 
=  O. 128353524 


=  0.4 


3  b , 


I  SEN 

TH. 

ACT. 


HERD  COEFF.  =  4.663401576 
i  E  G .  0  F  R  E  A  C  T 1 0  N  :=  3 .  0  51715  8  2  6 
DEG.  OF  REflCTIOH=  0„ 080315356 


'ELOCI 


i  '!' 


R I ANGLE  DATA 


',.,'  i  =  i  ■"' 
VA1  =  3 
VU1=  123! 

ALPHA  1= 
BETA  1  = 


10  ,,67141 
■  €  » 9  i  6  4  € 
i8.  22418 


Vfl2= 


7b. 40720 
6  O  ■  fc  9  2  4 13 


234,80211 
234.79915 
-1. 17934 

ALPHA  2=    -O 
BETH  2  =   -70 


6143 


POINT  #8 


PRESSURE  RAT  10 

REF  FLOW  RATE 

REF  RPM 

REF  ROTOR  MOMENT 

REF  HP 

EFF.  T-S 

ST ATOP  LOSS  THE OR. 
ST AT OR  LOSS  COEFF. 
ROTOR  LOSS  COEFF. 
ROTOR  LOSS  THEOR. 


=  2r9"4737£93 

=  1 .  0  2  8  0  3  0  0  3  O 

=  1  "'OSS,  73933 

=  137 . 6  8 9  0  4  4 

~  3  ?  .  2  b  ('  4  y  r'  9  6 

=  O. 1071O0795 

--■■  O.  139140155 

=  O.033159935 

=  0.414274733 


I  SEN.  HERD  COEFF.  =  4.169275243 
TH.  DEG.  OF  REACTION  =  0,033020949 
A  C  T .,  DEG,  0  F  R  E  A  C  T 1 0  N  =  O ,  O  9  6  8  5  9  9  7  2 

VELOCITY  TRIANGLE  DATA 


VI  =  120O.60325 
VA1=  298,71298 
VU1=  1  162:,  34940 


ALPHA  1= 
BETH  1  = 


"7  cr   i:~  n 
i  J  .  -J  .7 

c:;  q  ..  !=i  1  : 


'2  =  219.93553 
'A£=  2O3.69610 
'U2=    69, 43873 

ALPHA  2-         13 
BETA  2  =   --78 


3962  A 
v  1  ..;i-i!: 
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TABLE  III  TURBINE  TEST  RESULTS  (CONT.) 
POINT  #9 


P  R  E  S  S  U  R  E  R  A  T 10  =  2 .  5  1 9  4  8  8  0  3  4 

R E i::'  F L 0 W  P H  I E  =  1.8 2 8 099908 

REF  RPM  =  17176.36109 

REF  ROTOR  MOMENT  =  113.7867884 

REF  HP  =  31.00049275 

EFF.  T-S  =  0.743362397 

8  T  fl  T  0  R  L  0  S  S  "!"  H  E  0  R .  =  0 .  1  9  7 1 9  0  7  9  5 

S  T  R  T  0  P.  L  0  S  S  C  0  E  F  F  .  =  0.1 8  9  6  3  3  3  2 1 

R  0  T  0  R  L.  0  S  S  C  0  EFF,.  =  0.19129  2  8  5  7 

R 0  T 0 R  L  0  8 S  T H E 0  R .  =8.41 3 8 6  8 8 4  4 


I S E N .    H E  fl D    C 0 E  F  F .  =    3 .56765  6 1 2  2 

T H „    BEG.    0  F    R E  fl  C T 1 0  N    =   8.858 1 8 3 6  7  8 
fl C  T .    D E G .    0 F    R E fl C T 1 0 H -    8.18 8 3 8726 9 

VELOCITY  TRIANGLE  DATA 

VI  =  1103.46913  V2  =  249.75975 
Vfl  1  =  2  5  9 .  3  9  3  8  2  V  fl  2  =  1 8  3 .  8  4  8  2  4 
V U  1  =  1972.547  8 9    V U 2  =   169.959  8 8 

fl L  P H fl  1  =    7  6 . 4  8 4 1 8    fl  L  P  H  fl  2 =    4 £.6881 9 
BETH  1  =    57.94743    BETH  2  =   --69.39925 


POINT  #10 


P R E S  S  U R E  RATI  0  =  2.818415 7 9  8 

REF  FLOW  RATE  =  1.020999099 

REF  RPM  =  17153.62371 

REF  ROTOR  MOMENT  =  32.24633045 

REF  HP  =  22.39163977 

EFF.  T-S  =  9.683139969 

3 T  fl T Q R  L 0  S S  T H E 0  R .  =8.1 9  7 1 9  9  7  9  5 

8 T  fl  T  0 R  L  0  S 3  C 0  E F  F .  =  9.114113599 

R  0  T  0  R  L  0  8  8  C  0  EFF',.  =  9.4  8  9  2  3  8  8  9  8 

ROTOR  LOSS  THEOR.  =  9.351634739 

I  SEN,,  HERD  COEFF.     =  2.787522827 

TH.  DEC.  OF  REACTION  =  Q. 114732223 
ACT.  DEC.  OF  REflCTIOM=~l. 15697E-03 

VELOCITY  TRIANGLE  DATA 


VI    = 

9  8  7 .  9  9  9  8  0        V  2    =      3  4  5 .  2 1  7  8  5 

Vfll  = 

C  I  v  »  c  4  '-J  4  .1.            v  H ,.:.  -         i.  •_'  o  .  o  j  o  1  J 

VU1  = 

9  g  4  „  4  4  i  4  g         v  |  j  v  =      3  O  9 .  8  3  4  7  5 

ALPHA 

1=         77.78237         ALPHA    2=         S3.! 

■i  ~i  ?' ;:::: 

BETH 

■;        ._.               ETC       ,;:'  -|   .;•  ';::-•                nPTQ               ..-            ...  •"  i," 

'765 
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FIGURE    5       TURBINE  C 
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FIGURE  7    TURBINE  TEST  RIG  GEOMETRY  FOR 
TURBINE  CONFIGURATION  C 
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FIGURE  8   ILLUSTRATION  OF  SMOOTHING  USING 
FIRST  DEGREE  POLYNOMIAL 
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FIGURE  9   STATOR  LOSS  VS.  ISENTROPIC  HEAD 
COEFFICIENT,  BEFORE  SMOOTHING 
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FIGURE  10 
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STATOR  LOSS  VS.  ISENTROPIC  HEAD 
COEFFICIENT,  AFTER  SMOOTHING 
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FIGURE  11   ROTOR  LOSS  VS.  ISENTROPIC  HEAD 
COEFFICIENT,  BEFORE  SMOOTHING 
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FIGURE   12     ROTOR  LOSS   VS.    ISENTROPIC   HEAD 
COEFFICIENT,    AFTER  SMOOTHING 
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FIGURE  13   PREDICTED  PERFORMANCE  RANGE, 

RPM=15,000,  INPUT  ROTOR  L0SS=.25l4 
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FIGURE  14 


PREDICTED  PERFORMANCE  RANGE, 
RPM=18,000,  INPUT  ROTOR  L0SS=.25l4 
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FIGURE  15  PREDICTED  PERFORMANCE  RANGE, 

RPM=18,000,    INPUT  ROTOR   LOSS=.l 
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FIGURE  16  PREDICTED  REFERRED  HORSEPOWER  VS.  Kfel 
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FIGURE   17      PREDICTED  EFFICIENCY  VS.    Kfel 
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FIGURE  18   PREDICTED  STATOR  LOSS  VS.  K 
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FIGURE  19  PREDICTED  ROTOR  LOSS  VS.  K 
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FIGURE  20  REFERRED  HORSEPOWER  VS.  PRESSURE  RATIO 
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FIGURE    21   REFERRED  HORSEPOWER  VS.    PRESSURE  RATIO 
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MEASURED  REFERRED  HDR5EPDWER 
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FIGURE  22  VELOCITY  DIAGRAM  FOR  TU&BINE  TEST 
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FIGURE  23   VELOCITY  DIAGRAM  FOR  TURBINE  TEST 
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FIGURE  2k     VELOCITY  DIAGRAM  FOR  TURBINE  TEST 
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FIGURE  25  VELOCITY  DIAGRAM  FOR  TURBINE  TEST 
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FIGURE  26  VELOCITY  DIAGRAM  FOR  TURBINE  TEST 
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FIGURE    27      VELOCITY   DIAGRAM  FOR  TUftBINE   TEST 
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FIGURE    28      VELOCITY  DIAGRAM  FOR   TURBINE   TEST 


B 

X 

S 

r» 

— • 

r-x  li 

u 

SB 

•X. 

no. 

?i 

OS 

a. 

#1  6 

t- 

Q 

ZV.UJ 

—  bbq 

qi-Q 

OS 

a.  a.  a: 

a 

»- 


I 


mgnw 


r-mro 

wrnait/i 


ki- 
rn 


bin 


>-a 


in  — 

grni 


— a  1 1 

.  EC  — fM 

i  imiEo: 
3>>occca:aaj 


61 


FIGURE  29  VELOCITY  DIAGRAM  FOR  TURBINE  TEST 
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FIGURE  30  VELOCITY  DIAGRAM  FOR  TURBINE  TEST 
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FIGURE   31      VELOCITY  DIAGRAM  FOR  TURBINE   TEST 
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APPENDIX  A 
TURBINE  PERFORMANCE  ANALYSIS 

A-l   INTRODUCTION 

This  appendix  gives  a  detailed  description  of  the 
theory  and  analytical  technique  used  to  predict  the 
performance  of  a  turbine  in  the  Turbine  Test  Rig.   Variables 
used  in  this  appendix  are  defined  as  they  are  introduced. 
Also  given  are  the  details  of  the  computer  program  used 
to  implement  the  analytical  method.   The  variables  used 
in  the  computer  program  are  listed  and  defined  in  Table  A-I. 

A-2   ASSUMPTIONS 

The  following  assumptions  were  made: 

1.  The  stage  was  choked  in  the  stator. 

2.  In  solving  the  continuity  equation  between  the 
stator  and  rotor  the  solution  corresponding  to  the  lower 
static  pressure  was  taken. 

Assumption  (1)  was  made  on  the  basis  of  past  exper- 
ience with  the  test  rig.   Ref.  1  states  that  there  is 
supersonic  flow  in  the  stator  and  an  examination  of  the 
pressure  distribution  in  the  stator  passages  confirms  this 
statement.   Fig.  A-l  is  a  plot  of  the  pressure  distri- 
bution through  the  stator  by  pressure  tap  number  (See  Fig.  6). 
It  can  be  seen  in  this  plot  that  there  was  a  sharp  pressure 
rise  between  tap  #3  and  tap  #4,  which  indicates  the 
presence  of  a  shock  in  the  divergent  section  of  the  passage. 
Also,  the  level  of  static  pressure  at  the  throat  taps 
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in  comparison  to  the  supply  pressure,  was  observed  to  "be 
independent  of  downstream  pressure  for  all  operating  con- 
ditions. 

Assumption  (2)  limits  the  range  of  possible  solutions 
as  shown  in  Fig.  A-2.   This  figure  is  a  plot  of  flow 
coefficient  (0)   vs.  pressure  ratio  {?*/¥.)    for  various 
loss  coefficients  (z).   (The  flow  coefficient,  a  non-dimen- 
sional mass  flux,  is  introduced  in  Section  A-2.)   There 
are  two  possible  pressures  for  a  given  flow  coefficient 
and  given  loss  coefficient.   Solutions  were  limited  to 
the  left  side  of  the  line  connecting  the  points  of  maximum 
flow  function.   The  procedure  was  chosen  because  it  resulted 
in  the  "supersonic  root"  for  the  flow  from  the  stator. 
Also,  conditions  corresponding  to  choking  were  sought 
for  the  rotor,  and  the  locus  of  locally  sonic  conditions 
is  as  shown  in  Fig.  A-3.   It  should  be  noted  that  because 
of  the  definitions  of  the  flow  coefficient  and  loss  coef- 
ficient, the  "choking"  condition  at  which  the  downstream 
pressure  has  no  effect  upstream  of  the  plane  in  question, 
does  not  correspond  to  the  maximum  mass  flux  from  given 
stagnation  conditions. 

The  analysis  is  pseudo-1-dimensional.   Blockage  factors 
are  introduced  with  the  physical  cross  sectional  area  to 
account  for  non-uniform  flow  conditions.   Loss  coefficients 
are  defined  on  the  basis  of  kinetic  energies.   The  analysis, 
which  requires  only  mass  flow  continuity  through  the  stage, 
is  divided  into  three  stages;  the  stator,  the  interblade 
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region  and  the  rotor. 

While  the  analysis  itself  is  general,  solutions  can 
only  be  obtained  by  assigning  values  to  chosen  parameters. 
In  the  solutions  reported  here,-  the  flow  angle  at  the  stator 
exit  (<*j) ,  the  blockage  factors  for  the  stator  throat 
(kts)  and  rotor  exit  plane  (k,R),  and  the  rotor  passage 
loss  coefficient  (zR)  were  given  values.   The  value  ofCXj 
was  fixed  ato^  =75°   based  on  the  results  of  the  turbine 
tests  and  from  the  blading  geometry.   The  blockage  factors 
were  set  at  k*,s=.965  and  k,R=.85.   The  rotor  passage  loss 
coefficient  (zR)  was  an  input  variable  which  was  set  at 
.251^  (theoretical  value)  or  .1  (lowest  value  obtained  in 
test  results).   The  "rotor  loss  coefficient"  output  from 
the  program,  z-»_2»  includes  all  losses  from  the  stator  exit 
plane  to  the  rotor  exit  plane.   zi_?  corresP°nc*s  "k°  "the 
"rotor  loss  coefficient"  evaluated  from  test  rig  measure- 
ments* consequently,  zi_2^zR* 
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A-3   BASIC  RELATIONS 


The  sketch 
shows  a  general 
adiabatic  process 
with  entropy  increase 
(losses)  on  a  T-S 
diagram  on  which 
the  temperature 
scale  has  been  divided 
by  the  constant  stagnation  temperature  of  the  process.  A 
perfect  gas  is  assumed.  (0)  represents  the  initial  stag- 
nation state  and  (  )  represents  the  final  state. 

The  non-dimensional  velocity,  X,  is  defined  as 

X  =-JL/f  where  Vt  =V§c"T^is  the  "total"  or  "limiting" 

Vto  P 

velocity,  and  the  loss  coefficient,  z,  is  defined  as  shown. 

Note  that  the  stagnation  velocity  is  constant  throughout 

the  process. 

The    "flow  coefficient",   0,    is   defined  here   as 


0m 


J^l 


A(l) 


where  w  is   the   flow  rate,       Pftg is   the    density  at  the   initial 
stagnation  temperature   and  pressure,   A   is   the    flow  area  and 
kfe  is   the    blockage   factor,      The   flow  coefficient  defined 
in  this  way  is   a  non-dimensional  flow  rate   per   unit   area,    or 
mass   flux,    referred  to   initial  stagnation  conditions. 

The   flow  rate    is   given  by 

WmfAkyV 
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so   that  Eq.   A(l)    can  be   written 

which,    in  terms   of  the   non-dimensional  velocity  defined 
above,    becomes 

The    loss   coefficient  is   defined,    as   shown  in  the    above 
sketch,    as 

z  =     T    ~7*'*   =    /  -      -*£-  Aj>3) 

Tto  -Tis  Xt$. 

Rearranging  Eq.   A(3)i 


*-*)•-£--     -^ 


x%        I-  p- 

and  using  the  isentropic  relationship  between  pressure  and 
temperatures 

Using  Eq.  A(*0,  Eq.  A(2)  becomes  a  single  equation 
for  0   as  a  function  of  pressure  ratio  P/P+0»   In  analysing 
the  turbine,  generally  the  flow  rate  is  known  so  that  the 
flow  coefficient  can  be  calculated  using  Eq.  A(l).   On 
specifying  a  value  for  the  loss  coefficient, ' (z) ,  the  pres- 
sure ratio  can  be  calculated  from  Eq>  A(2)  using  Eq.  A(^). 
An  iterative  technique  is  used  (Newton's  Method)  starting 
with  an  initial  estimate  of  the  pressure  ratio.   Fig.  A(2) 
shows  0   as  a  function  of  P/P+Q  for  various  values  of  zt 
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in  graphical  form.   The  solution  to  the  left  of  the  maxima 
is  obtained  by  beginning  the  iteration  with  a  small  value 

of  VPto- 

The  analysis  of  the  turbine  begins  by  specifying  that 
the  flow  rate  is  set  by  the  choking  of  the  stator  nozzles. 
The  pressures  at  successive  stations  are  then  calculated 
in  turn  from  the  flow  coefficient  as  described  above. 

A-4  STATOR  EXIT  CONDITIONS 

The  flow  coefficient  at  the  throat  of  the  stator  is  the 
value  of  the  flow  coefficient  corresponding  to  locally 
sonic  conditions  at  the  minimum  area.   In  the  present  work, 
the  maximum  flow  coefficient  at  a  loss  coefficient  equal  to 
.05  (defined  as  0*)   was  taken  at  the  stator  throat, 
throughout  the  calculations. 

Since 


*  =  r  v  TTt  A(5) 


where  A  and  k,  are  the  area  and  blockage  factor  at  the 
throat  station,  the  flow  rate  being  constant  through  the 
machine  requires  that 


^  =  fu   Vt.   4  k9s  0*   =  constant 


at  stations  ahead  of  the  rotor,  where  T.   is  constant  J  this 

to 

implies   that 

P6ff  fts  I*ls0*  =   constant 
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The   value   of  0*   can  be   obtained  analytically  when  the 
loss   coefficient   is   known.      The   maximum  occurs  where 
~j" .=   0,    which  leads   to   the   solution   (denoting  values   at  the 
stator  throat  by  an  asterisk), 

(X*)2  =   B  -/b2  -c'  A(6) 

where   B  =    YV  *  '*      »    C  =    (l-z^-JrM 
to      1  j  +  p 


and  P*/P+rt   ■  (-rPz    )  A(7) 


whereD=^/^£l 


Then  0*   is   given  by  Eq.   A(2). 

With  0*   known,    continuity  of   flow  rate   requires   that 
at   the    stator  exit, 

d   =  0*  Ain A(8)      I 

where  subcript  1  represents  the  stator  exit  plane  (station  1) 
and  ©<, is  the  flow  angle  at  thestator  exit. 

Knowing  0*  ,  P-j/P+o  can  be  found  by  iteration  as 
described  in  section  A-3.   Solutions  on  the  left  side  of 
Fig.  A-2  were  selected  by  beginning  the  iteration  with 
P. /P.   close  to  zero. 

After  obtaining  thevalue  for  P-j/P+q  then  X.  can  be 
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calculated  using 


\  -h-^ 


h 


P,  \  » 


fro 


A(9) 


Then,  Pi/P-j^  can  be  found  fr 

Vpti  -  (1-xi2)'"' 


om 


A(10) 


and 


Ptl/ptc  =  <VPto>  <*W 


A(ll) 


A-5   INTERBLADE  SPACE 

The  interblade 
space,  from  Station 
(1)  to  Station  (i) , 
is  shown  in  the  sketch. 

Conditions  at  (i) 
are  calculated  in  a  manner  -   //  /   / 
similar  to  those  at  Station  (1). 

First,  by  continuity, 


Akk. 


£ 


1      ^i    h^Cost;         Pt, 


A(12) 


where  &1   is  the  flow  angle  at  the  rotor  entrance. 
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The  value  of  c<i    is  obtained  from  the  condition  that  ang- 
ular momentum  is  conserved  in  the  interblade  space.   Then 


X 


R 


Ul     R 


1-Xj  sino<,  A(13) 


where  R.=mean  radius   at  Station  1,    and  R.=mean  radius   at 

Station  i. 

Then 


'■  -   .«„-!   Xui 


°<l  =  sin"1  -JiL  k(lk) 

A  • 

i 


and,    also 


\~(Vr-*)A.  ^ 


where   N  is   the   RPM,    and  the   tangential  velocity   (V   . ) 


has   been  non-dimensionalized  as 


to 

Knowing  0.  ,  a  value  for  P;/p+i  can  be  found  if  a 
value  of  the  interblade  loss  coefficient  (z.)  is  assumed. 
Then 


VPto  ■  <VPtl>   (ptl/pto'  A<17) 
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and 


X. 

1 


A(18) 


With  the  absolute  flow  properties  established  at 
the  rotor  entrance,  the  relative  flow  conditions  can  be 
calculated.   First,  from  the  geometry  of  the  velocity 
diagram,  the  relative  flow  angle  (^ )  is  given  by 


&  - 


=   tan 


-1 


Xc  sia/<x; :-U: 


At  COfi  <*t 


A(19) 


where 


U. 


.  Ui 


to 


A(20) 


is  the  non-dimensional  rotor  speed  at  radius  R. .   Then  the 
non-dimensional  relative  velocity  (X  . )  is,  from  the  vel- 

Wl 

ocity  diagram,  given  by 


X   =  cos  **«- 
wi    cos  Qi 


A(21) 


and  the  temperature  T. ,  by 


T.        l 
to 


A(22) 


7^ 


The  equivalent  temperature  (Ref.  1),  is  given  by 


-El  =  1  -  x2  +  x2.  +  u2 

T  l     wi     1 

ito 


and  therefore 


PEi   _      PEi   Pi 


to 


P. 

l 


to 


Ri 


-  1 


TMh^-'  pi 


to      1 


to 


A(23) 


A(2^) 


where 

Tt-..    =  equivalent  temperature    into  rotor    (See  Fig.   A-?) 

£1 

PEi  =  equivalent  pressure  into  the  rotor  (See  Fig.  A-?) 

A-6   ROTOR  CONDITIONS 

Using  the  values  calculated  thus  far  it  is  now  possible 
to  calculate  the  conditions  in  the  rotor.   Continuity 
requires  that 


Pe. 


A(25) 


W 


where  a  subscript  R  denotes  the  exit  plane  in  the  rotor 
blading. 

It  should  now  be  noted  that  the  process  in  the  rotor 
frame  from  the  equivalent  conditions  at  (i)  to  the  exit 
of  the  rotor  is  entirely  similar  to  the  basic  process 
described  in  Section  A-3»   However,  the  equivalent  temper- 
ature  takes  the  place  of  stagnation  temperature  and 
relative  velocity  replaces  velocity  in  the  given  equations. 
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The  pressure   ratio  which  satisfies   the    flow  coefficient 

given  by  Eq.    A(25)f    for   an  assumed  rotor  passage   loss 

coefficient    (z,-,)    is  Po/P^. .      This   is   obtained  as   described 
R  2     Ei 

in  Section  A->.      The   corresponding  non-dimensional  velocity 
is   now 


^-.M-h) 


-  Cl. 


id. 


£c 


A(26) 


where 


W. 


w2 


PTiT 


A(2?) 


Hence 


X...^  =  Y 


•Ei 


w2  w2./  T 


to 


A(28) 


From  the   velocity  diagram, 


X2  =    (ir2   *   Xw2  sin?J2+    <Xw2   C0S^)2 


A(29) 


where  Q%    =  71°  (from  the  geometry  of  the  rotor),  and 


2   ui  RL 
The  temperature  at  the  rotor  exit  is  given  by 


to 


•) 


E: 


<1-«k)  [T-r  zr 


?to. 


so  that 


■^2    h-    +   x2 
Tto  "   Tto    X2 


A(30) 
A(3D 
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Next  the  loss  coefficient  that  includes  all  of  the 
losses  in  the  stage  from  Station  1  to  Station  2  can  be 
calculated.   This  coefficient  includes  all  of  the  losses  in 
the  inter-blade  region  as  well  as  those  in  the  rotor 
passage.   It  is  called  z.,  2  here  but,  in  fact  corresponds 
to  the  measured  rotor  loss  coefficient  in  the  TTR  since  all 
losses  aft  of  the  stator  in  the  TTR  are  included  in  the 
"rotor  loss  coefficient." 


,         VT2 
z<    «=1- 


1-2  TE-T2, 

IS 


71 


s  1 Zi« Ttf 

Tie            Tz;t 
=    1 ?**■ 2i- ^  t  A(32) 


A    Pi 


i   rt,\  y 


XzL 


A-7  STAGE  PERFORMANCE  PARAMETERS 

In  the  notation  used  for  the  turbine  test  rig  in 
Ref.  1  and  Ref.  2,  the  following  equations  determine  the 
turbine  performance  parameters  from  quantities  calculated 
in  the  above  steps  1 

V    Tt2  \ 

^w-'tol1  -T£"J  A(33) 


AT.      =  T. 
is    to 


2  \  r 


A(3f) 


^T-S  "  S^  A(35) 


^Tis 
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p     vJto  toAo  bo; 


H.P.    =  w   (,2*K>2)  -?HATW 


M  =  H.P.    (550) 


360 


«r= 

p 

6    - 

T 

T 
iref 

V>     = 

H.P. 

*    -      H.P. 

M*   = 

M 

N*    = 

N 

^ 


A(36) 

A(37) 
A(38) 

A(39) 
A(^0) 
A(M) 

A(42) 

A(43) 
A(4^) 
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A -8   THE  COMPUTER  PROGRAM 
A-8.1  Overview 

There  are  nine  degrees  of  freedom  in  the  com- 
puter program: 

1.  C^i   ,  the  flow  angle  out  of  the  s  tat  or 

2.  k,.  ,  the  stator  exit  blockage  factor 

3.  ^bi»  ^e  interstage  blockage  factor 

k.      k,R,  the  rotor  passage  blockage  factor 

5.  z  ,  the  stator  loss  coefficient 

6.  z.,  the  interstage  loss  coefficient 

7.  zR,  the  rotor  passage  loss  coefficient 

8.  N,  turbine  RPM 

9.  T.  ,  total  temperature  upstream  of  the  stator 
o<,  is  input  on  the  basis  of  turbine  test  results. 

At  high  pressure  ratios  the  value  of  &,   was  measured  to 
be  75°-76°.   Therefore,  the  value  of  ©(,  was  set  and  held 
at  75° '•   The  values  of  RPM  and  T.   are  those  conditions 
for  which  a  solution  is  being  sought. 

All  the  other  parameters,  the  rotor  passage  loss 
coefficient  (zR)  and  blockage  factor  (k,R)  were  given 
chosen  values.   This  choice  was  made  because  the  incidence 
losses  were  included  in  the  interblade  calculation.   It 
was  thought  that  the  parameters  describing  the  flow  inside 
the  rotor  passages  could  be  held  fixed  as  the  speed 
was  allowed  to  vary. 

Consequently,  there  are  four  degrees  of  freedom  with 
which  the  program  works.   The  goal  for  the  program  is  to 
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find  those  solutions  which  result  in  flow  coefficients  in 
the  rotor  that  are  less  than,  or  equal  to,  the  "choking" 
value.   This  is  done  in  the  program  by  iteration  on  the 
remaining  four  degrees  of  freedom.   Fig.  A-5  is  a  block 
diagram  showing  schematically  the  iteration  process.   The 
order  of  iteration  is  (from  the  most  frequent  to  the  least 
frequent),  z.,  z  ,  k,  .  ,  and  k,  .  .   The  loss  coefficients  for 
both  the  stator  and  interblade  always  begin  at  zero  and  they 
are  incremented  as  the  program  seeks  solutions.   A  graphical 
representation  of  a  typical  iteration  is  shown  in  Fig.  A-4. 
It  can  be  seen  that  the  program  can  find  a  solution  for 
the  condition  where  the  loss  coefficient  is  equal  to  zero 
(denoted  as  (1)  on  Fig.  A-4) .   On  the  next  iteration  the 
program  will  set  the  loss  coefficient  equal  to  .1  and  it 
can  be  seen  that  a  solution  still  exists  at  point  (2).   On 
the  next  iteration,  however,  it  can  be  seen  that  a  solution 
does   not  exist  for  a  loss  coefficient  equal  to  .2.   In  this 
case  the  program  will  go  back  and  add  a  smaller  increment 
to  the  previous  loss  coefficient  (.1). 

Each  time  an  added  increment  results  in  a  loss  coef- 
ficient greater  than  the  maximum  allowable  the  program 
subtracts  that  increment  and  adds  a  smaller  one.   In  this 
way  the  program  eventually  finds  the  maximum  loss  coef- 
ficient which  will  yield  a  solution.   At  this  point, the  given 
flow  coefficient  is  at  the  maximum  point  for  the  calculated 
loss  coefficient  line.   At  this  point  "choking",  as  defined 
here,  has  occurred. 
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It  should  be  pointed  out  that  this  "choking"  con- 
dition (the  point  of  maximum  flow  coefficient  for  any- 
given  loss  coefficient)  is  not  the  point  corresponding 
to  locally  sonic  conditions.   In  Fig.  A-3  the  locus  of 
points  of  maximum  flow  coefficient  and  the  locus  of  points 
of  local  sonde  flow  are  both  shown.   The  point  of  locally 
sonic  flow  always  falls  to  the  left  of  the  "choking"  point 
(as  defined  in  this  work),  and  the  flow  coefficient  for 
sonic  flow  is  always  less  than  the  flow  coefficient  at 
the  defined  choking  point.   The  reason  is  that  the  point 
of  maximum  mass  flux  (maximum  flow  coefficient)  is  for 
givsn  upstream  stagnation  conditions  (see  next  section). 
Because  of  the  way  the  flow  coefficient  is  defined,  the 
maximum  for  a  given  loss  coefficient  will  not  correspond 
to  local  sonic  conditions  erxcept  where  the  loss  coefficient 
is  equal  to  zero. 

The  difference  between  the  "choking"  flow  coefficient 
corresponding  to  the  maximum  and  the  flow  coefficient  at 
the  local  sonic  condition  is  about  10$  at  a  loss  coefficient 
of  .3  (the  highest  rotor  loss  coefficient  predicted  by 
the  program).   The  difference  is  about  k%   at  a  loss  coef- 
ficient equal  to  .25  (the  highest  stator  loss  predicted), 
and  insignificant  (less  than  1%)    at  a  loss  coefficient  of 
about  .15. 
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A-8.2  Description  of  the  Program 

The  following  discussion  refers  to  Fig.  A-6, 
which  is  a  detailed  flow  chart  of  the  program,  to 
Table  A-I ,  which  is  a  list  of  program  variables  and  their 
definitions,  and  to  the  program  listing  given  in  Section  A-10. 

The  first  step  in  the  program  is  to  input  the  variables. 
c(,  is  input  in  line  280  (it  is  also  converted  to  radians 
in  that  line).   kfe1  is  input  in  line  286  and  k,  .  is  input 
in  line  28^-.   All  other  variables  are  input  as  requested 
by  the  program  (lines  20-110).   Note  that  P?  and  P.   are 
input  also.   They  are  not  involved  in  the  calculations  and 
were  originally  input  as  the  conditions  for  which  a 
particular  s/olution  was  sought. 

The  next  step  is  to  go  to  a  subroutine  to  calculate 
the  value  of  0*    (line  170).  0*    is  defined  as  the  value 
of  the  flow  coefficient  at  the  throat  of  the  stator  which 
sets  the  values  of  the  flow  coefficients  at  all  downstream 
stations.  0*   is  generated  in  the  subroutine  by  the  method 
covered  in  Section  A-^  assuming  a  loss  coefficient  of  .05 
to  the  throat  of  the  stator  (line  160). 

Next  the  program  calculates  the  conditions  in  the 
stator  exit  plane.   This  is  done  by  first  calculating  the 
value  for  0.    (line  320).   As  explained  earlier,  the  value 
for  0*    determines  the  upper  limit  of  the  stator  loss 
coefficient. 

Using  the  method  described  in  Section  A-3»  "the  program 
calculates  the  value  for  P-/P+  .   Note  that  if  the  stator 
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loss  coefficient  is  large  enough,  so  that  a  solution  does 

not  exist,  this  will  be  detected  in  the  decision,  "P./P.  <0 

1   to 

or>l"  (lines  350  and  370).   If  "yes"  then  the  stator  loss 

(z_)  is  too  large  and  the  following  steps  are  takeni 
s 

1.  The  last  increment  added  to  the  stator  loss  is 
subtracted  (line  4-20).  (Note,  the  initial  stator  loss 
increment  is  .1 . ) 

2.  The  existing  stator  loss  increment  is  multiplied 
by  some  factor  less  than  1  (line  bj>0) . 

3.  The  new,  smaller,  increment  is  then  added  back 
to  the  stator  loss  coefficient  (line  1570). 

If  stator  loss  increment  is  sufficiently  small 
(line  1575)  then  the  upper  limit  of  the  loss  coefficient 
has  been  reached  as  discussed  earlier.   At  this  point  no 
further  solutions  are  attainable  and  k,  .  is  incremented 
(line  2290).   Note  that  the  upper  limit  of  stator  loss 
will  not  normally  be  attained  until  several  iterations 
have  been  made  through  the  entire  program. 

Next  the  interblade  calculations  are  made.   The  sol- 
ution for  P.1/P.1.  found  for  the  stator  exit  plane  is  used 
to  make  the  calculations  for  0.    (line  770).   Again,  0. 
determines  the  upper  limit  on  the  interblade  loss.   In 
the  interblade  region,  however,  there  is  another  criterion 
that  must  be  met.   The  value  for  ^;  generated  within  the 
program  must  be  equal  to  £  (g  =69°  and  is  fixed  by  the 
geometry).   This  criterion  is  met  by  iterating  the  inter- 
blade loss  coefficient  as  follows  1   On  every  pass  through 


83 


this  section  the  interblade  loss  is  initially  set  to  zero 
and  the  interblade  loss  increment  is  initially  set  to  .1 
(lines  720,  7^0).   The  value  for  P-/Pt«  is  calculated  using 
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same  procedure  as  described  earlier  for  P-i/P+0«   Again, 
the  value  of  P. /P..  is  checked  on  each  iteration  to  ensure 
that  it  falls  between  0  and  1  (line  780).   If  it  does  not 
then  the  interblade  loss  is  reduced  in  a  method  analagous 
to  that  covered  for  the  stator  loss. 

After  generating  a  value  for  Pj/P-m  "the  program  cal- 
culates §i   (line  1208).   When  the  value  for  Q    is  generated 
it  is  compared  to  ^g  and  three  results  are  possible. 
If  p«  is  equal  to  Q    ±.05  then  converganne  has  occurred  and 
the  program  will  continue.   If  ^  is  less  than  Q     (line  1250) 
then  the  value  for  interblade  loss  coefficient  is  increased 
by  the  current  value  of  interblade  increment  and  the  program 
goes  back  to  the  start  of  the  interblade  calculations  and 
begins  anew.   If  ^.  is  greater  than  @ff   (line  12^0)  the 
value  of  interblade  loss  is  decreased  in  the  same  way 
as  was  described  for  the  stator  loss. 

Two  results  are  possible  in  the  interblade  regiont 
Either  the  program  finds  the  value  of  the  interblade  loss 
that  gives  flow  angle  convergence  ( ^;  =  £  ) »  or  the  inter- 
blade loss  increment  becomes  sufficiently  small  (10~  )  to 
trigger  termination. 

Note  that  the  interblade  loss  increment  is  allowed  to 
become  very  small  (10~  )  (line  850).   This  is  because  the 
above  mentioned  angle  convergence  requirement  is  very 
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sensitive.   Note,  also,  that  when  the  interblade  loss 
increment  falls  to  10~  and  the  iteration  process  in  the 
interblade  region  is  terminated,  the  stator  loss  coefficient 
is  incremented  (line  850).   Stator  loss  is  incremented 
"because  a  higher  value  of  stator  loss  may  result  in  con- 
ditions that  allow  convergence  in  the  interblade  region. 
If  not,  the  stator  loss  calculations  eventually  reach  the 
termination  criterion  mentioned  earlier. 

Once  angle  convergence  has  been  attained  in  the  inter- 
blade region  then  the  value  of  the  rotor  flow  coefficient 
can  be  calculated  (line  1^70).   If  the  program  gets  to  this 
point  for  any  given  set  of  conditions  then  a  solution 
exists  and  the  only  determination  to  be  made  is  whether 
the  resulting  rotor  flow  coefficient  is  greater  than  the 
"choking"  value  (line  1530).   If  Myes"  then  it  is  dis- 
allowed, since  a  flow  coefficient  greater  than  the  "choking" 
value  cannot  occur  physically.   If  the  calculated  coef- 
ficient is  less  than,  or  equal  to  the  "choking"  value  then 
the  results  are  printed,  the  stator  loss  is  incremented,  and 
the  calculations  start  again  at  the  stator  exit  plane. 

The  choking  value  of  the  flow  coefficient  for  the 
rotor  (#*R)  is  determined  in  exactly  the  same  way  as  the 
choking  value  for  the  stator  (0*),   with  the  assumed  value 
of  rotor  loss  used  as  the  loss  coefficient.   (At  present 
the  value  for  0t   is  calculated  separately  and  inserted  in 
line  121.   It  is  suggested  that  a  "GOSUB  2140"  be  inserted 
at  this  point.   Then  the  input  value  for  rotor  loss 
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coefficient  can  be  changed  easily  and  the  "choking"  value 
for  the  rotor  will  be  calculated  in  the  program.)   If 
the  resulting  flow  coefficient  is  greater  than  "choking", 
the  program  increments  stator  loss  coefficient  and  returns 
to  the  stator  plane  calculations  to  begin  anew. 

Therefore,  in  all  cases  the  program  will  seek  the 
highest  value  of  stator  loss  coefficient  for  which  solutions 
exist.  At  some  point,  the  rotor  flow  coefficient  normally 
attains  a  value  less  than  the  "choking"  value  and  all 
solutions  for  higher  stator  loss  coefficients  will  be 
acceptable  solutions.   Or,  if  the  "choking"  condition  is 
never  reached,  the  program  will  attain  the  highest  possible 
stator  loss  in  attempting  to  achieve  solutions  less  than 
"choking" . 

When  the  highest  value  of  the  stator  loss  is  reached, 
there  are  two  possible  avenues  in  the  program.   If  the 
program  was  generating  solutions  less  than  "choking"  at 
the  time  of  termination  in  the  stator  calculations,  then 
k,  .  will  be  incremented,  stator  loss  will  be  set  back  to 

Dl 

zero,  and  the  process  begun  anew  (line  2280).   If  the 

solutions  being  generated  at  the  time  of  termination  were 

greater  than  "choking",  then  the  upper  value  of  k,  .  has 

been  reached  and  any  increase  in  k,  .  will  yield  no  further 

solutions.   In  this  case  k,  „  is  incremented,  k,  .  is  set 

bl  bi 

to  the  pre -determined  lower  limit,  and  stator  loss  is 
reset  to  zero  (line  2330) • 

The  program  continues  in  the  above  manner  until  the 
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pre -determined  upper  limit  of  k,..  is  reached  (line  2335)  i 
at  which  point  the  program  is  terminated. 

Since  the  21-MX  computer  has  only  6  digits  of  accuracy, 
round-off  error  has  been  a  problem.   Where  Newton's  Method 
is  used  to  solve  for  pressure  ratio  the  computer  lacks 
sufficient  accuracy  to  converge  to  the  desired  accuracy 
under  certain  conditions.   When  this  happens,  the  program 
will  stay  inside  the  iteration  and  never  converge.   This 
problem  has  been  solved  by  putting  a  counter  inside  each 
iteration  loop  (lines  5^0  and  980).   More  than  60  iterations 
is  treated  as  a  no n- convergent  condition,  and  the  approp- 
riate loss  coefficient  is  reduced  in  the  manner  covered 
previously. 

A-9   OPERATING  PROCEDURE 

The  procedures  for  operating  the  Hewlett-Packard 
21-MX  computer  will  not  be  covered  in  this  paper  since 
the  applicable  manuals  are  available  at  the  laboratory. 
The  computer  must  be  on,  with  the  RTE-B  operating  system 
"READY". 

1.  Load  the  paper  tape  program  labeled  "TTR  11". 

2.  Edit  into  the  program  the  minimum  values  of  k,. 
(line  28*0  and  k,.  line  (286).  Edit  the  desired  value  of 
o(t  into  line  280.   If  it  is  desired  to  get  a  print-out 

of  every  solution  then  either  remove  lines  172*4-  and 
1726  or  place  "REM"  in  front  of  them.   If  these  two  lines 
are  left  in  the  program  then  a  solution  will  be  printed 
out  only  if  the  calculated  pressure  ratio  (Pp/^t  ^  ^s  ^-ess 

87 


than,  or  greater  than,  the  previously  calculated  lowest  or 
highest  pressure  ratio  respectively,  for  the  value  of  k,  1 
then  being  used. 

3.   Type  "RUN",  "RETURN"  on  the  keyboard  and  the 
program  will  begin  execution  by  asking  for  inputs: 

"INPUT  STATOR  LOSS  COEFFICIENT"  The  normal  input  is 
zero.  However,  if  some  specific  case  is  required  then  any 
desired  loss  coefficient  can  be  input. 

"INPUT  PTO"  Input  the  upstream  total  pressure. 

"INPUT  P2"  Input  the  hood  pressure. 

"INPUT  RPM"   Input  the  RPM  desired. 

"INPUT  TTO"  Input  the  total  temperature  into  the 
stage. 
This  is  all  that  is  required  to  operate  the  program. 
Depending  on  the  range  of  k,  .  requested  the  program  will 
take  from  30  minutes  to  36  hours  to  run  on  the  21-MX. 

It  is  important  to  realize  that  the  above  procedure 
will  produce  solutions  that  give  values  of  the  rotor 
flow  coefficient  less  than  "choking".   If  the  turbine 
performance  at  the  "choking"  condition  is  desired  then 
after  the  computer  run  is  complete  a  further  step  must  be 
taken.   It  is  necessary  to  take  the  values  of  k,  .  and  k,  . 
for  which  a  solution  is  desired  and  force  the  choking 
condition.   This  is  done  as  follows* 

1.   Scan  the  output  results  at  the  desired  k,.  and 

k,  .  and  locate  the  point  that  has  a  value  of  rotor  flow 

Dl  r 

coefficient  closest  to  the  choking  value  and  a  non-zero 
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stator  loss  coefficient. 

2.  Re-start  the  program,  this  time  putting  in  a  value 
of  the  stator  loss  coefficient  slightly  less  than  the  value 
printed  on  the  output. 

As  stator  loss  is  manually  decreased,  the  rotor  flow 
coefficient  will  increase  towards  the  choking  value.   When 
the  rotor  flow  coefficient  is  within  the  limits  specified 
in  line  2370  the  words  "CONVERGENT  CONDITIONS"  will  appear 
on  the  teletype  output  preceding  the  printed  results. 

If  the  stator  loss  coefficient  is  decreased  too  much, 
the  rotor  flow  coefficient  will  increase  to  a  value  above 
choking.   In  this  case  no  ouput  will  be  printed  on  the  tele- 
type.  However,  the  flow  coefficient  will  appear  on  the 
video  display.   In  this  case  increase  the  stator  loss 
slightly. 

It  has  been  found  in  this  work  that  the  extreme  values 
of  referred  horsepower,  stator  loss,  rotor  loss,  and 
efficiency  occur  at  the  lowest  and  highest  values  of  k,  . 
for  each  value  of  k,  1 .   With  this  in  mind  it  is  possible  to 
find  the  range  of  predicted  values  at  choking  by  forcing 
to  choking  only  the  lowest  value  of  k,  .  (with  non-zero 
stator  loss)  and  the  highest  value  of  k,.  for  each  k,  . . 
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A-10    PROGRAM  LISTING 
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TABLE  A-I 
DEFINITION  OF  VARIABLES 

AO  4A*)   Area  at  the  Stator  Throat 

Al  (A.)   Area  at  the  Stator  Exit  Plane 

A8  Decision  Variable 

A9  Decision  Variable 

CI  Decision  Variable 

D2      dUL, 

D3  (ATW) 

E  Decision  Variable 

EO  ^t-S^    Total-to-Static  Efficiency 

FO  (0*)      Flow  Coefficient  at  Stator  Throat 

Fl  (jZL)      Flow  Coefficient  at  Stator  Exit  Plane 

F2  (0R)     Flow  Coefficient  Through  Rotor  Passage 

F3  (A)      Flow  Coefficient  Through  Interblade   Area 

F^  Input  Value   of  0     Based  on  Assumed  Rotor  Loss 

F8  Decision  Variable 

F9  Decision  Variable 

GO  (Gamma  -   1)/Gamma 

HO  Horsepower 

HI  Referred  Horsepower 

H8  {§)   Rho 

H9  Reynolds   Number 
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j  or)    pi 

Kl  k?s  Blockage  Factor  at  Stator  Throat 

K^  k,  .  Interblade  Blockage  Factor 

K5  k,  R  Blockage  Factor  in  the  Rotor 

K6  k,  ,  Blockage  Factor  at  the  Stator  Exit  Plane 

Ml  Rotor  Torque 

M2  Referred  Rotor  Torque 

N  RPM 

NO  Referred  RPM 

01  (#,)   Flow  Angle   at  Stator  Exit  Plane 

03  (^t)   Flow  Angle   at  Rotor  Entrance   Plane 

0^-  (Ak,.  )      Increment  for  k,  . 

P  Pti/Pto 

P0  Pto 

P1  VPto 

P2  VPto 

p3  Pl/ptl 

P*  Ptl/pto 

p5  pi/ptl 

P6  VPto 

p7  PEi/pto 

P8  Decision  Variable 

P9  P2 

R  (fto)      Total  Density  at  Stagnation  Conditions 

R7  Relaxation  Parameter 

R8  Relaxation  Parameter 


R9     Predicted  P2/pt0 
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S7-S9  Used  in  GOSUB   7000  to  Find  P0/P-,. 

C.  HI 

TO  TtQ 

T1  Tt2/Tto 

T2  T/T^ 

T3  T./Tt0 

T*  TEi/Tto 

T5  TtAt0 

T6  WTto 

Ul  U.  Dimensionless  Velocity  (See  Appendix  A) 

U2  U"2  Dimensionless  Velocity  (See  Appendix  A) 

U3  U.  Dimensionless  Velocity  (See  Appendix  A) 

VO  V.   Total  Velocity 

WO  w  Flow  Rate 

Wl  w*  Referred  Flow  Rate 

W9  Decision  Variable 

XO  X  .  Non-Dimensional  Relative  Velocity 

XI  X.   Non-Dimensional  Velocity 

X2  X2  Non-Dime nsional  Velocity 

X3  X.   Non-Dime nsional  Velocity 

X^  X  .  Non-Dime nsional  Velocity 

X5  X  .  Non-Dime nsional  Relative  Velocity 
wi 

X6  X  p  Non-Dimensional  Relative  Velocity 

X8  Decision  Variable 

Y2  Non-Dimensional  Velocity 

Y6  Increment  for  k,  . 

Y7  Increment  for  k,. 

Y8-Y9  Counters 
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ZO  z  Stator  Loss  Coefficient 

Z3  z.  Interblade  Loss  Coefficient 

Z?  zi-2  Predicted  Rotor  Loss  Coefficient 

Z8  Increment  for  Stator  Loss  Coefficient 

Z9  Increment  for  Interblade  Loss  Coefficient 
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FIGURE  A-l   PRESSURE  DISTRIBUTION  THROUGH  THE  STATOR 
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FIGURE   A-2      ILLUSTRATION  OF  SOLUTION  AREA 
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FIGURE  A-3   MAXIMUM  0   AND  SONIC  LINE 
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FIGURE  A-4-  TYPICAL  SOLUTION  PROCESS 
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FIG.  A-5  PROGRAM  SCHEMATIC 
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FIGURE  A-?  THERMODYNAMIC  PROCESS  OF  FLUID 

IN  AN  AXIAL  TURBINE  STAGE 


APPENDIX  B 
TURBINE  TEST  RIG  (TTR)  DATA  REDUCTION  AND  PROCESSING 

B-l   INTRODUCTION 

This  appendix  describes  only  those  changes  that  have 
been  made  to  the  data  reduction  procedure  given  in  Ref.  2. 
The  program  numbers  and  the  functions  of  those  programs 
given  in  Ref .  2  have  not  changed.   The  channel  and  port 
assignments  for  data  collection  and  storage  have  not 
changed.   Variables  also  have  not  changed,  although 
additional  variables  were  defined  in  the  course  of  mod- 
ifying the  data  reduction  process.   Those  variables  that 
now  exist  in  addition  to  those  given  in  Ref.  2  are  given 
in  Table  B-I . 

Modifications  made  during  the  course  of  this  work  were* 

1.  Raw  data  can  now  be  read  directly  from  mass  memory. 
Previously  it  was  necessary  to  read  the  paper  tape  for  a 
given  point  in  order  to  reduce  it. 

2.  Provisions  have  been  written  into  the  programs 
for  smoothing  the  reduced  data. 

3.  All  of  the  parameters  used  in  calculating  the 
theoretical  loss  coefficients  that  were  previously 
entered  from  charts  are  now  in  polynomial  form  in  the 
program. 

4.  Storage  of  raw  data  is  now  a  separate  program. 
Previously  it  was  necessary  to  run  the  reduction  sequence 
to  store  data. 
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B-2   DESCRIPTION  OF  PROGRAMS 

The  data  reduction  is  divided  into  ten  separate 
programs  with  an  additional  program  to  store  the  raw  data. 
Fig.  B-l  shows  the  contents  of  each  program  and  also  shows 
the  sequence  in  which  the  programs  are  chained.  Following 
is  a  description  of  the  reduction  sequence  referred  to 
Fig.  B-l  which  will  be  followed  by  a  description  of  how  to 
run  the  reduction  program. 

TTR  in  Fig.  B-l  is  used  only  to  store  the  raw  data. 
TTR  does  not  chain  to  any  other  program. 

The  data  reduction  process  begins  in  TTR1.   In  order 

to  smooth  the  reduced  data  as  discussed  in  Section  III  it 

is  necessary  to  have  n  values  for  (P.  ~PhuOAPti  ~Phub^  • 

Therefore,  the  first  step  in  the  data  reduction  process  is 

to  go  from  TTR1  to  the  point  in  TTR2  where  the  values  of 

P-t  P+'-rj*  and  Puu-ut  have  been  calculated.   On  completion, 

the  program  returns  to  TTR1  and  repeats  the  process  n 

times  until  the  n  values  of  (P. -P.  ,  )/(P+.  -P.  .  )  have 

1  hub    tip  hub 

been  calculated.   The  program  then  proceeds  to  TTR9  where 
a  polynomial  curve  fit  for  the  n  values  of 
(P1-PhubHPti  -Phub)  is  generated.   Then  the  program 
returns  the  polynomial  coefficients  to  TTR1B  where  the 
reduction  process  begins  at  the  first  point  in  the  run. 
However,  when  the  point  in  TTR2  is  reached  where  P1/P+  is 
calculated,  it  is  calculated  using  the  polynomial  as 
described  in  Section  III. 

The  rest  of  the  reduction  process  in  Fig.  B-l  is  the 
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same  as  described  in  Ref.  2  with  the  exception  that  no 
further  keyboard  inputs  are  needed  until  the  raw  and 
reduced  data  is  tabulated. 

B-3  RAW  DATA  STORAGE 

The  procedure  for  storing  raw  data  on  the  mass  memory 
is  outlined  in  Ref.  2,  Appendix  B.  Sec.  B-3»  lines  10?19. 
Note  that  line  10  should  read,  "GET  ,TTR,H.   Although  it  is 
now  a  separate  program  the  procedures  for  storing  data 
have  not  changed. 

B-k     PROCEDURE  FOR  DATA  REDUCTION  PROGRAM 

1.  Key  in  GET  "TTR1B" 

2.  Press  "RUN  EXECUTE". 

3.  The  following  check  list  will  be  printed  on  the 
HP  9830  printer » 

"PRIOR  TO  RUNNING  THIS  SEQUENCE  ENSURE  FOLLOWING i 
"THAT  TTR2  LINE  590  HAS  PROPER  FACTOR  IN  IT; 

1.02  FOR  HOODED 

1.01   FOR  UNHOODED 

"IF   IT   IS   NECESSARY  TO   INPUT  ALPHA    1   THEN  CHANGE   TTR2 

"LINE   960  TO    'INPUT  A31 

"IF   BLOCKAGE  FACTOR   OTHER  THAN   1,    CHANGE   TTR2-1060 

"TO    'INPUT   X7',    AND  PUT  SEMICOLON  AFTER   TTR2-1050 

"ENSURE   THAT   TTR1 ,    TTR1B  HAVE   THE   PROPER   FILENAME   IN   W*l 

"IF   IT   ISN'T   DESIRED  TO  STORE   REDUCED  DATA   CHANGE 

TTR6-580  TO    'G1=0' 

"ENSURE  TTR6-610,  63 0  HAVE  PROPER  FILENAME  FOR  REDUCED 

DATA 
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"ENSURE  TTR7-220  HAS  PROPER  FILENAME  FOR  RAW  DATA 
"ENSURE  TTR8-220  HAS  PROPER  FILENAME  FOR  REDUCED  DATA 
"IF  IT  IS  DESIRED  TO  RUN  A  SINGLE  POINT,  WITHOUT 
SMOOTHING,  INSERT  TTR2-1127  'GOTO  1200'.   OTHERWISE 
OMIT  TTR2-1127 
The  above  checklist  will  prepare  all  ten  programs.   The 
last  item  is  a  provision  for  elimination  of  smoothing  if 
it  is  desired  to  run  a  single  point.   If  smoothing  is  not 
desired  then  each  point  must  be  run  separately  with 
TTR2-1127  inserted. 

b.      Press  CONTINUE  EXECUTE 

5.  The  display  will  read  LOWEST,  HIGHEST  RECORD  # 
THIS  RUN.   Input  the  lowest  record  number  and  the 

highest  record  number  from  which  it  is  desired  to  read 
raw  data.   Note,  the  record  number  on  which  the  reduced 
data  will  be  stored  will  be  the  same  number  as  the  raw 
data  record  from  which  the  raw  data  was  read.   Make  sure 
the  filename  designated  for  the  storage  of  the  reduced 
data  can  accept  data  on  those  record  numbers  without 
writing  over  previous  information. 

6.  Next  the  display  will  read  PRINT  OUT  RESULTS? 
YES=1,  N0=0.   If  a  1  is  selected  here  then  at  the  conclusion 
of  data  reduction,  for  all  the  points  in  the  run,  the 
program  will  begin  to  print  all  results.   If  0  is  selected 
then  all  reduced  data  will  be  stored  on  the  mass  memory. 

7.  It  will  take  about  15  minutes  for  the  program  to 
read  the  raw  data  and  make  the  calculations  necessary  for 
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smoothing  the  data.   When  this  is  oompleted  the  program  will 
start  reducing  each  point  and  printing  the  results. 

8.  At  the  completion  of  data  reduction  for  all 
points  the  option  to  tabulate  the  data  is  available. 

ENTER  RECORD  # 'S :  LOWEST,  HIGHEST  will  appear  on  the  display. 
If  it  is  desired  to  tabulate  the  raw  data  then  enter  the  same 
record  numbers  that  were  entered  to  initiate  the  program. 
The  raw  data  will  then  be  tabulated. 

9.  Next  the  display  will  read  ENTER  LOWEST,  HIGHEST 
RECORD  NUMBER.   If  it  is  desired  to  tabulate  the  reduced 
data  then  enter  the  appropriate  record  numbers.   The 
reduced  data  will  be  tabulated. 
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TABLE  B-I 
VARIABLES  ADDED  TO  DATA  REDUCTION  PROGRAM 

A2  Lowest  Record  Number 

A5  Counter 

A6  Decision  Variable 

A7  Highest  Record  Number 

B  (Array)  Polynomial  Coefficients  From  TTR9 

F8  Decision  Variable 

Q  (Array)  Values  of  Isentropic  Head  Coefficient 

T^+  Polynomial  Approximation  of  Reynolds  Number 

Z  Decision  Variable 
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